Collective electromagnetic relaxation in crystals of molecular magnets 
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We study the magnetization reversal and electromagnetic radiation due to collective Landau- 
Zener relaxation in a crystal of molecular magnets. Analytical and numerical solutions for the time 
dependence of the relaxation process are obtained. The power of the radiation and the total emitted 
energy are computed as functions of the crystal parameters and the field sweep rate. 
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I. INTRODUCTION 

Paramagnetic crystals of high-spin molecular magnets, 
like Mni2, Fes an d others, exhibit unusual magnetic 
properties related to the macroscopic time of the tran- 
sition between spin-up and spin-down states of individ- 
ual magnetic molecules jlj. The latter is due to the high 
magnetic anisotropy and a large value of spin, S 3> 1. 
For, e.g., a biaxial molecule (Fes of S = 10), in the mag- 
netic field, H, parallel to the anisotropy axis Z, the spin 
Hamiltonian is 



H 



-DS 2 + AS 2 - gf i B H z S z 



(1) 



where g is the gyromagnetic factor, /j, B is the Bohr mag- 
neton, and D > A > 0. For small A, the approximate 
energy states of H are the eigenstates of S z : S z \m) — 
m\m). At H z = kD/g^ B , with k = 0, ±1,±2, ..., the 
levels m < and m! satisfying m + rn! = —k come 
to resonance. For the even S, the tunnel splitting of 
the resonant levels, A m , appears (for even k) in the 
[{m! — m)/2]-th order of the perturbation theory on A: 
A (m= _ s) ex (A/D)( m '- m )/ 2 . At, e.g., k = (see Fig.QJ, 
A( m= _ S ) oc (A/D)° and, thus, at S = 10, the probability 
of the transition between spin-up and spin-down states is 
low. Consequently, at low temperature, the crystal can 
be prepared in a state with inverse population of the spin 
energy levels, e.g., magnetized against the direction of the 
magnetic field. This allows one to observe, in a macro- 
scopic experiment, such quantum effects as resonant spin 
tunneling J3, Q > s P m Berry phase Q , crossover between 
quantum tunneling and thermal activation y| la, , and 
quantum selection rules in the absorption of electromag- 
netic radiation [8|. 

Recently, it has been suggested 0, EH E3 that a crys- 
tal of molecular nanomagnets can be a source of co- 
herent electromagnetic radiation in the millimeter wave- 
length range, highly desirable for applications [T^- Some 
experimental evidence of this effect has been obtained 
M [l^, Ull • The effect is related to Dicke superradiance 
[l5j . Normally, atoms or molecules of a gas, initially pre- 
pared in the excited energy state, decay independently 
by spontaneous emission of light. The power of the radi- 
ation obeys the law P ex N exp(— t/r) where N is the 
total number of atoms and r is the lifetime of the excited 
state. Dicke argued that N atoms confined within a vol- 




FIG. 1: Approximate energy levels of a spin- 10 molecule in 
a zero magnetic field. The tunnel splitting of the degenerate 
levels is not shown. Arrows show the relaxation path from 
m — — 10 to m = 10 through thermally assisted quantum 
tunneling. 



ume of size d , which is small compared to the wavelength 
of the radiation A, cannot radiate independently from 
each other. At d < X a spontaneous phase locking of 
the atomic dipoles takes place, that results in the coher- 
ent radiation burst of power Psr oc N 2 , emitted within 
a time of order tsr ~ t/N. This phenomenon, called su- 
pcrfluorescence, has been widely observed in gases. It can 
occur in any system of identical quantum objects if the 
system is not very large compared to the wavelength of 
the radiation [l|| . For crystals of molecular magnets this 
is true for both, the transitions between the tunnel-split 
levels, Fig. and the transitions between the adjacent 
\m) levels, Fig.Q 

In a typical experiment, one magnetizes the crystal and 
then sweeps the field in the opposite direction. In this pa- 
per we will be concerned with the situation when the elec- 
tromagnetic transitions occur between tunnel-split levels, 
Fig. These can be, e.g., transitions between m = —10 
and m = 10 levels shown in Fig. ^ The electromag- 
netic relaxation of that kind corresponds to the total 
magnetization reversal accompanied by the broad band 
superradiance. It is described by a rigorous model [2( 




FIG. 2: A pair of tunnel-split levels vs. energy bias W. The 
Landau-Zener (LZ) transition is followed by the emission of 
the coherent light via superradiance. 



which is reviewed in Section II. In essence, if one ne- 
glects the electromagnetic radiation, the crossing of the 
(to, to') resonance by the magnetic field sweep, Fig. is 
described by the Landau-Zener theory [l7| . When the 
coupling between the spins and the electromagnetic radi- 
ation is taken into account, the magnetic state resulting 
from the LZ transition relaxes towards the lowest energy 
state via superradiance. The rate of the superradiant 
decay, as well as the time dependence of the relaxation, 
are sensitive to the parameters of the crystal and to the 
shape of the magnetic field pulse. Our goal is to com- 
pute the time dependence of the radiation power and the 
total radiated energy as functions of the field-sweep rate, 
the tunnel splitting, and the size of the crystal. This is 
done in Section III by analytical and numerical meth- 
ods. Practical implications of our findings are discussed 
in Section IV. 



II. COLLECTIVE LANDAU-ZENER 
RELAXATION 

Consider a crystal of N magnetic molecules occupying 
an to magnetic state that is close to the resonance with 
the to' state, e.g. to = —S, ml = S in Fig.^ We shall as- 
sume that the molecules weakly interact with each other 
through the electromagnetic field. As has been shown 
in Ref. |9j , the quantum magnetic relaxation of such a 
crystal satisfies the Landau-Lifshitz equation: 



h = 7[n x H cff ] - cry [nx[nxH cff ]] 



(2) 



Here n is a unit vector of the pseudospin describing the 
two-state system, such that n z — —1 corresponds to all 
molecules in the TO-state, while n z = 1 corresponds to all 
molecules in the m'-state, 7 = g[i B /H is the gyromagnetic 
ratio, the effective magnetic field is given by 




a = 0, e = 0.1 



FIG. 3: Time dependence of the magnetization reversal for 
two values of e due to pure Landau-Zener relaxation of indi- 
vidual molecules (a = 0) and due to collective relaxation via 
superradiance (a = 0.01). 



with W{t) — \gn B {m' — m )H z being the energy sweep 
rate (see Fig.|2j|, and a <C 1 is a dimensionless effective 
damping coefficient, 



^N(m' - m) 2 g 2 a ( 



m e C' 



(4) 



with a = e 2 /he ss 1/137 being the fine structure con- 
stant. Note that a is independent of the magnetic field. 
The first term in Eq. @ gives dissipationless Landau- 
Zener transitions when the field is swept through the res- 
onance such that W = W(t) satisfies W(±oo) — ±00, 
and the initial condition is n(— 00) = — e z . Indeed, at 
a = the Schrodinger equation for a two-level system 
is equivalent to the equation for a precessing spin. The 
probability p(t) for the molecule to stay in the initial 
state is given by 



p(t) = [1 - n,(t)]/2 



(5) 



For W(t) — vt, one obtains the famous Landau-Zener 
result [17]: p(oo) = p LZ = cxp(-e), where 



irA 2 
2hv 



(6) 



The Landau-Zener effect corresponds to only partial 
magnetization reversal, 



LZ 



(00) = 1 — 2exp(— e) 



(7) 



g^ B H c s = Ae^ + We z 



(3) 



see Fig. at a = 0. The value of n z z (oo) is close to 
-1 at e < 1, that is for the fast field sweep. In this 
case most of the molecules, after crossing the resonance, 
remain in the initial m-state by passing from the lower 
to the upper branch in Fig. [5J On the contrary, for a 
slow sweep, that is when e> 1, most of the molecules 
follow the lower branch in Fig. [21 and the final state of 
the crystal is exponentially close to n z z (oo) = 1. 

The second term in Eq. J5J) describes collective mag- 
netic relaxation via Dicke superradiance. Due to this 
term the magnetization of the entire crystal at long times 



reverses completely to n z (oo) = 1, as is shown in Fig. 
for a finite a. The collective relaxation due to superradi- 
ance is significant for e < 1, that is, when n z z (oo) is not 
very close to 1. Thus, the observation of the superradi- 
ance requires a fast field sweep. 



III. RADIATION POWER 

The power of the superradiance described by Eq. (J2J) 
can be obtained from the classical formula for the mag- 
netic dipole radiation |9j]: 



where 



P(t) = [2/(3c")]mi(t) 



m z (t) = -N(m - m')gn B n z (t) 



(8) 



(9) 



Close to the (m, m!) resonance, nearly any field sweep of 
practical interest is linear in time, W — vt. It is conve- 
nient to use dimensionless variables: 



t' = 



tA ,, , s vt hvt' 7rf' , N 



In terms of these variables Eq. (J2J and Eq. (JHJ) become 



<7n 
and 



7 = [n x (e x + W'{t')e z )] - a[nx [nx (e x + W'{t')i 



P = aNh~ 1 A x 



j2 \ 2 

an. 



dt' 



(11) 
(12) 



The total emitted energy, E — J dtP(t), is given by 

E = aNAE', (13) 

where we have introduced dimensionless 



E' = / dt' 



J2 x 2 

drn z 



dt 1 



(14) 



A. Analytical 

We shall start by developing an analytical approxima- 
tion for the practical case of e < 1 and a<l. The time 
interval of interest is the one past the Landau-Zener tran- 
sition: W 1 ^> 1. In this case, retaining the leading terms 
in Eq. (JTJJ, we get 



(15) 
(16) 
(17) 



dn x 


= W'n y 


dn y 
IF 


= -W'n 


dn z 
dt' 


= -% + 



-sf = -n v + aW'(t')(l-ni) 
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FIG. 4: Approximate analytical solution for n z (t') averaged 
over oscillations, Eq. (121)1 . for two values of e and a — 0.01 
(solid line). The numerical solution of Eq. II U is shown by 
the dash line. 



These equations show that n x and n y oscillate rapidly in 
time, while n z , in accordance with Fig. |3 has a slowly 
varying average. Averaging Eq. I|17|l over the period of 
oscillations of n v , one obtains: 



dn z 
IF 



= aW'{t')(l-nl) 



(18) 



Eq. H18fl describes the superradiant stage of the evolu- 
tion of n z . Therefore, it must be solved with the initial 
condition n z = n^ z at t — 0. At small e, Eq. Q gives 
for that initial condition: 

n z (0) = -l + 2e. (19) 

The corresponding solution of Eq. (|18|l reads 14] 



/ n , I a7rt 1 , 1 

n z (t = tanh In — 

y ' \ 4e 2 e 



(20) 



It is shown by the solid line in Fig.0] For e = 0.1 Eq. fTI^I 
is, clearly, a good approximation to the full solution aver- 
aged over oscillations. As e increases, some discrepancy 
is observed. One can improve the analytical approxi- 
mation by writing n z — h z + Sn z and solving Eq. I|17|) 
through iterations, but this, at the end, will require a 
numerical integration, so that the improvement obtained 
by this method does not give much advantage over the 
direct numerical solution of Eq. (|ll|l . 

The approximate solution of Eqs. (|15|) and i|ltj|) . satis- 
fying n 2 = 1, is given by 



h> = \A - n 2 z sin ( — 



-VT 




(21) 
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FIG. 5: Time dependence of the reduced radiation power, 
P' = (d 2 n z /dt' 2 ) 2 at e = 0.3 and a = 0.05. Solid line repre- 
sents numerical results. Dash line corresponds to Eq. I|22> at 
<A„ = 2.576. 



FIG. 6: The e dependence of the total emitted energy at 
two values of a. Points represent numerical results. Solid line 
corresponds to Eq. 1231 . 



where <fi Q is a phase which we are not attempting to com- 
pute analytically. One can see from Eqs. (|2U|I and i|21|) 
that n z , indeed, changes slowly with time, compared to 
the oscillations of n x (t') and n y (t'), because of the con- 
dition a<l. 

Let us now turn to the analytical approximation for the 
power and the total radiated energy. It is easy to see from 
Eqs. (|17|I - 1|21JI that the main contribution to d 2 n z /dt' 
is determined by the rapidly oscillating n^-term in the 
right-hand-side of Eq. (fl 71) : 



d n, 



dt 



l2 



tt£_\ sin[(7rf'74e) + (j> ] 
2e ) cosh[(a7rf' 2 /4e) + ime] 



(22) 



Substituting this expression into Eq. (|14|l and replacing 
the rapidly oscillating sin 2 (irt' /4e) under the integral by 
1/2, one finally obtains: 



E' = 



V2 a 3/2y o cosh 2 (a .2 + l lne ) 



(23) 



Equations JEJ, CSJ, GU, (E3, and (J23J) give the de- 
pendence of the radiation power and the total emitted 
energy on the field sweep rate and on the parameters of 
the crystal. 



B. Numerical 
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We shall now compute P(t) and E by numerical inte- 
gration of Eq. (|ll|l , and compare them with our analyti- 
cal findings. 

The time dependence of the reduced power, P' = 
(d 2 n z /dt' ) 2 , is shown in Fig. The comparison with 
Eq. lt22l) is performed by fitting the value of <j) a until where 
the match with the numerical solution of Eq. (|11|) for 
(d n z /dt' ) 2 is obtained. Even for e as large as 0.3 the 



FIG. 7: The a dependence of E at two values of e. Points 
represent numerical results. Solid line corresponds to Eq. 
P3)i. 



agreement of the numerical results with the analytical 
formula is rather good. Note that the oscillation of the 
power in time is a quantum effect related to the oscilla- 
tion of n z . 

Fig. HO shows the dependence of the total emitted en- 
ergy on the parameter e, that is, on the inverse field sweep 
rate. Fig.[7]shows the dependence of E' on the parameter 
a. 

The question of significant importance for experiment 
is the spectral composition of the radiation. The total 
emitted energy can be presented as 



E= diuI(Lu) 



I(co) = haNI'(uj') 



(24) 



(25) 



IV. DISCUSSION 
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FIG. 8: Spectral function I'{ui') for three values of alpha at 
e = 0.1. 
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The formulas and the numerical results obtained above 
are valid if the energy distance between the resonant 
levels, W, is small compared to the distance between 
the adjacent m-levels. For numerical estimates, we shall 
stick to the (—5, 5)-resonance. The conclusions of this 
section, however, will apply to other resonances as well. 
For the model illustrated by Eq. (|TJ and Fig. ^ the dis- 
tance between the m — —5 level and the m = —5 + 1 
level is (25 - l)D sw 25L>. The validity condition we 
are looking for is then W(t) < 2SD. For W = vt 
one should verify this condition at W max — vt max , 
where t max — ht' max /A is the time when the superra- 
diance drops exponentially due to the hyperbolic cosine 
in Eq. I|22|) . According to this equation and Eq. I|10(l . 
t'ma X ~ \/ e / a an d Wma X ~ A/ ' \[ea. Substituting this 
into Wmax <C 2SD and using Eq. Q for a, one obtains 
the validity condition in the form of the lower bound on 
the total number of molecules: 



1 f m e c 
ea V U a 



(28) 



0.05 0.1 0.15 0.2 0.25 0.3 



FIG. 9: Dependence of / of Eq. J23 on e. 



is the spectral power. Here I'(u>') is a dimensionless func- 
tion of the dimensionless frequency, lu' = fko/A. It must 
be computed via the Fourier transform of d 2 n z /dt' : 



Z7T 



dt'e v 



dt' 



(26) 



This function is shown in Fig. [S] The peak of the power 
occurs at 



huj 



ma X — i — J \£) 



(27) 



This scaling of hw max on a follows from Eq. (|llfl . The 
function /(e), computed numerically, is shown in Fig. 



where U a — DS 2 is the energy barrier between m = ±5 
states due to magnetic anisotropy. Eq. I|28|l shows that 
a high magnetic anisotropy and a not very small e are 
needed if the size of the system is to remain within rea- 
sonable limits. The optimal would be e ~ 1 since, ac- 
cording to Eq. (JZJ), e < 1 (that is, a sufficiently high 
field-sweep rate) is needed to create an inverse popula- 
tion of spin levels. For Mni2 and Fes, the anisotropy bar- 
rier is of order 60K and 30K respectively, and the lower 
bound on N, according to Eq. I)28|) . must be between 
10 18 and 10 19 molecules. With account of the unit cell 
volume (3.7 nm 3 and 2.0 nm 3 for Mni2 and Fes, respec- 
tively) this translates into a volume of order or greater 
than 1mm 3 . Remarkably, this agrees with the reported 
lower bound on the volume of the crystal (or crystal as- 
sembly) that shows evidence of electro mag netic radiation 
during magnetization reversal [T(J, ES 03 • 

We shall now estimate the total emitted energy and 
the power of the radiation. According to Eq. 1)23(1 , E' ~ 
e" 1 /^" 3 / 2 . This gives for E of Eq. (JT3J: E - NA/y/ea. 
With the help of Eq. {2J one obtains: 



E 



-^N^mrC 2 



(29) 



For the purpose of the order-of-magnitude estimate we 
have dropped the factor gS^fao of order unity. No- 
tice that the total emitted energy is proportional to 
the square root of the crystal volume. At e ~ 1 and 
N ~ 10 18 , Eq. J2HJ) provides E - 0.1 mJ. 

According to Eqs. |Q| and l(2"2"|) (see also Fig. [SJ the 
power of the radiation oscillates in time. In most cases, 
observation of these oscillations must be impeded by the 
finite time resolution of the measuring equipment, so that 
only the envelope of the curve shown in Fig. [S] will be 
observed. The peak power can be estimated as P max ~ 



obtains 



NA 2 /he. Substituting here e of Eq. ©, one 



Nv ~ 5M- 



,dH 
dt 



(30) 



where we have introduced SM = g/i B (m' — m)N, the 
change in the total magnetic moment of the crystal due 
to collective electromagnetic relaxation. Note that the 
relations E ex y/~N and P ma x oc N are specific to the 
radiation problem we have studied. For an assembly of 
a few mm-size crystals Eq. <|3"U)l gives P max ~ 10 uW at 
a typical laboratory field sweep rate of 0.01 T/s [!(], Il3l 
and P max - 1 W for a fast field pulse, dH/dt ~ 10 3 T/s 
|l4| . Note, however, that in the case of an ultrafast sweep 
the condition e ~ 1 can be satisfied only by a large tun- 
nel splitting A, making the preparation of the initially 
magnetized state less simple than in the case of small A. 
During the adiabatic sweep, the frequency of the ra- 
diation is determined by the distance between the spin 
levels, hw — VA 2 + W 2 . For the most of the re- 
laxation process, W 7$> A, and, thus, u> = W{t)/K. 
The peak of the spectral power, Fig. [8J corresponds to 
hiOmax ~ W (t max ) ~ A/yfea. Up to a factor of order 
unity, that depends logarithmically on e, this coincides 
with Eq. i|27|) . The logarithmic difference of /(e) in Fig. 
El from l/y/e is due to hie in Eq. (J22J . With the help of 
Eq. J3J, we obtain that by order of magnitude 



m e c 



U), 



(31) 



where we again omitted the factor gS \/ao of order unity. 
For e ~ 1 and N ~ 10 18 , required to produce significant 
radiation (see above), this frequency is in the terahertz 
range. 

For the radiation to be coherent, the inhomogeneous 
broadening of e must be small throughout the crystal. 



This translates into narrow distribution of the tunnel 
splitting and narrow distribution of the magnetic field 
felt by the spins. Both conditions must be satisfied in 
Fes because they are also the necessary conditions for the 
Berry phase effect observed in that system Q. In Mni2 
the situation is less clear due to solvent disorder, large 
hyperfine interactions, dislocations, etc., which result in 
a distribution of A EH IS H^ EL [^ The suitability 
of Mni2 for the study of superradiance depends on 
whether the distribution of A is continuous or consists 
of a finite number of narrow lines due to, e.g., finite 
number of nuclear spin states, presence of isomers in the 
structure of the molecule |20||. etc. Systems with narrow 
distribution of A probably exist among hundreds of new 
molecular magnets synthesized in recent years. In such 
systems, narrow distribution of the magnetic field should 
be achieved automatically when the spins of the initially 
saturated sample rotate coherently due to superradiance. 

V. CONCLUSIONS 

We have studied magnetic relaxation via coherent elec- 
tromagnetic radiation, produced by the magnetic field 
sweep in a crystal of molecular magnets on crossing the 
tunneling resonance. We find that the effect exists start- 
ing roughly with crystals (or crystal assembly) of mil- 
limeter size. The radiation is broadband with the cutoff 
in the terahertz range. The power of the radiation is 
proportional to the field-sweep rate and ranges from mi- 
crowatts to watts for the existing sweep rates. 
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